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Abstract 

The theory of symmetries of systems of coupled, ordinary differential equations 
(ODE's) is used to develop a concise algorithm for cartographing the space of 
solutions to vacuum Bianchi Einstein's Field Equations (EFE). The symmetries 
used are the well known automorphisms of the Lie algebra for the corresponding 
isometry group of each Bianchi Type, as well as the scaling and the time reparam- 
eterization symmetry. Application of the method to Type III results in: a) the 
recovery of all known solutions without prior assumption of any extra symmetry, 
b) the enclosure of the entire unknown part of the solution space into a single, 
second order ODE in terms of one dependent variable and c) a partial solution to 
this ODE. It is also worth-mentioning the fact that the solution space is seen to 
be naturally partitioned into three distinct, disconnected pieces: one consisting of 
the known Siklos (pp-wave) solution, another occupied by the Type III member of 
the known Ellis-MacCallum family and the third described by the aforementioned 
ODE in which an one parameter subfamily of the known Kinnersley geometries 
resides. Lastly, preliminary results reported show that the unknown part of the 
solution space for other Bianchi Types is described by a strikingly similar ODE, 
pointing to a natural operational unification as far as the problem of solving the 
cosmological EFE's is concerned. 
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1 Introduction 



Since the early times of cosmology, Automorphisms have been identified as possible key 
elements for a unified treatment of spatially homogeneous Bianchi Geometries pQ. Har- 
vey has found the automorphisms of all 3-dimensional Lie Algebras while the corre- 
sponding results for the 4-dimensional Lie Algebras have been reported in Jantzen's 
tangent space approach sees the automorphic matrices as the means for achieving a 
convenient parametrization of a full scale-factor matrix in terms of a, desired, diagonal 
matrix Samuel and Ashtekar were the first to look upon Automorphisms from a 
space viewpoint (HJ. The notion of Time-Dependent Automorphism Inducing Diffeomor- 
phisms (A.I.D.'s), i.e. coordinate transformations mixing space and time in the new 
spatial coordinates and inducing automorphic motions on the scale-factor matrix, the 
lapse and the shift has been developed in 0. 

In this communication we revisit the problem of solving the EFE's for vacuum Bianchi 
Geometries. We begin with a full metric, i.e. we make no assumption for the lapse 
function iV 2 , the shift vector N a and the spatial metric y a p. Then we use the Time- 
Dependent A.I.D.'s to put the shift vector to zero. At this point the idea is to exploit, 
in a systematic way, the remaining symmetries of the field equations -sometimes called 
"rigid" [T.- to transform them to the most simple form possible, without loss of gener- 
ality. These are the well known symmetries following from the constant Automorphism 
group within each Bianchi Type, as well as the scaling of the metric by a constant and 
the time reparameterization symmetry (see e.g. jHj). Applying this analysis to Bianchi 
Type III Vacuum Cosmology we produce an exhaustive cartography of the entire space 
of its solutions. 

The paper is organized as follows: in section 2, we present our method. In section 
3, after a brief description of Bianchi Type III Cosmology we apply the method. We 
thus recover all known solutions, describe the unknown part of the solution space with a 
single, second order ODE in terms of one dependent variable and present a new solution. 
A brief preview of the corresponding results obtained for other Bianchi Types is also 
included. Finally some discussion and concluding remarks are given in section 4. 



As it is well known, for spatially homogeneous spacetimes with a simply transitive action 
of the corresponding isometry group jTOJ , [Hj, the line element, assumes the form 



2 The Method 



ds 2 = (N a N a - N 2 ) dt 2 + 2N a a?dx i dt + -f a ^a^dx i dx j 



(2.1) 



where the 1-forms erf, are defined from: 



da a = A cr 7 afj 

Then the field equations are (e.g. 0): 



(2.2) 
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(2.3) 



2 



E a = K»CL - K?C £ = 



(2.4) 



E a p = K aP + N (2K T a K Tf3 - KK a p) + 2N P {K av C v Pp + - NR a p = (2.5) 

where 

K aP = ~ (j aP + 2 lav C v p p N» + 2^ u C» p No) (2.6) 
is the extrinsic curvature and 

R Q/3 = CZrC^wpxr?" + 2Q A Ci + 2C^ xl , u r X + 

(2.7) 

zchflbrtairr + 2C* K c^ x r u 

the Ricci tensor of the hyper-surface. 

In jH| particular spacetime coordinate transformations have been found, which re- 
veal as symmetries of (|2.Hjl . ()2.4j) . (|2.5jl the following transformations of the dependent 
variables N, N a , 7^ : 

N = N, N a = A£ (iV p + 7 p. P"), %u = ^^lap (2.8) 
where the matrix A and the triplet P a must satisfy: 

K C M = C ^ A ?K (2-9) 



2P"C"A^ = (2.10) 

For all Bianchi Types, this system of equations admits solutions which contain three 
arbitrary functions of time plus several constants depending on the Automorphism group 
of each type. The three functions of time, are distributed among A and P (which also 
contains derivatives of these functions). So one can use this freedom either to simplify 
the form of the scale factor matrix or to set the shift vector to zero. The second action 
can always be taken, since, for every Bianchi type, all three functions appear in P a . 

In this work we adopt the latter point of view. When the shift has been set to zero, 
there is still a remaining "gauge" freedom consisting of all constant A~ (Automorphism 
group matrices). Indeed the system (|2.9J) . (j2.10|) accepts the solution A^ = constant, 
P a = 0. The generators of the corresponding motions, induced in the space of dependent 
variables spanned by 's (the lapse is given in terms of 7 Q/ g, 7 Q/ g by algebraically solving 
the quadratic constraint equation ) 7^ = A^ A^ y a p are [TT] : 

x m = y ma% ,^- (2.11) 

with A satisfying: 

^(I)p C/3-y = tfl)/3 C PJ + Vh C PP' ( 2 - 12 ) 
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Now, these generators define a Lie algebra and each one of them induces, through 
its integral curves, a transformation on the configuration space spanned by the 7 a( g's. If 
a generator is brought to its normal form (e.g. t£-), then the Einstein equations, writ- 
ten in terms of the new dependent variables, will not explicitly involve zi. They thus 
become a first order system in the function z% ^2J- If the above Lie algebra happens 
to be abelian, then all generators can be brought, to their normal form simultaneously. 
If this is not the case, we can diagonalize in one step the generators corresponding to 
any eventual abelian subgroup . The rest of the generators (not brought in their normal 
form) continue to define a symmetry of the reduced system of EFE's if the algebra of 
the Xm's is solvable [13]. One can thus repeat the previous step, by choosing one of 
these remaining generators. This choice will of course depend upon the simplifications 
brought to the system at the previous level. Finally if the algebra does not contain any 
abelian subgroup, one can always choose one of the generators, bring it to its normal 
form, reduce the system and search for its symmetries (if there are any). Lastly, two 
further symmetries of (|2.Hjl . (|2.4jl . (|2.5jl are also present and can be used in conjunction 
with the constant automorphisms: The time reparameterization t — > f(t) + a, owing to 
the non-explicit appearance of time in these equations, and the scaling by a constant 
7a/3 - ► as can be straightforwardly verified. Their corresponding generators are: 

(2.13) 

(2.14) 

well as with the Im's, as it can be 



y - 1 9 

Y 2 = 7«/3 ' 



<97a/3 

These generators commute among themselves, as 
easily checked. 



3 Application to Bianchi Type III 

We are now going to apply the Method, previously discussed, to the case of Bianchi 
Type III. For this type the structures constants are ^1] 

n\ _ _r<i — i 

°13 — °3i — 1 (3 1) 

= for all other values of a(3j 

Using these values in the defining relation ()2.2j) of the 1-forms <rf we obtain 

/0 e~ x 0\ 

of = 1 (3.2) 




The corresponding vector fields ^ ( satisfying = cTg£ 7 ) with respect to which 

the Lie Derivative of the above 1-forms is zero are: 

Zi = d y & = d z & = d x + yd v (3.3) 
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The Time Depended A.I.D.'s are described by 

A 



"73 



e- 2P W x(ty 

C 22 C 23 

1 



P ( 



x(t)P(t) + -x(t),P 2 (t),P(t) 



(3.4) 



(3.5) 



where P(t),x(t) and P 2 (t) are arbitrary functions of time. As we have already remarked 
the three arbitrary functions appear in P a and thus can be used to set the shift vector 
to zero. 

The remaining symmetry of the EFE's is, consequently, described by the constant 
matrix: 



Sl 



M 



S4 

e S2 s 3 
1 



(3.6) 



where the parametrization has been chosen so that the matrix becomes identity for the 
zero value of all parameters. 

Thus the induced transformation on the scale factor matrix is 7^ = M£Mg'j lll/ , 
which explicitly reads: 

7n = e 2si 711 
712 = e Sl+S2 712 



713 = e Sl (s 3 711 + s 4 712 + 713) 



722 = e 232 722 



723 = e S2 (a 3 712 + s 4 722 + 723) 



(3.7) 



, 733 = S3 2 711 + 2 S3 (S 4 712 + 713) + S 4 2 722 + 2 S 4 723 + 733 

The previous equations, define a group of transformations G r of dimension r = 
dim(Aut(1 1 1)) = 4. The four generators of the group, can be evaluated from the relation: 



X A 







9s A J s=0 dlafj 

where A = {1, 2, 3, 4}. Applying this definition to ()3.7|) we have the generators: 

d d d 

Ai = 2711- h 712^ h 713 q — 

0711 d7i2 0713 



(31 



(3.9) 



X3 = 712 7T 1 - 722 



7ir 



'0713 

d 



1X2- 



0723 
d 



27 



d 



23 



+ 2 7 



0733 



0713 0723 0733 

The algebra g r that corresponds to the group G r has the following table of commu- 
tators: 

[X 1; X 2 ] = 0, [Xl,X 3 ] = 0, [Xl,X 4 ] = x 4 , , . 

[x 2 ,x 3 ] = x 3 , [x 2 ,x 4 ] = o, [x 3 ,x 4 ] = o 

As it is evident from the above commutators ()3.13j) the group is non-abelian, so we 
cannot diagonalize at the same time all the generators. However, if we calculate the 
derived algebra of g r , we have 

g r , = {[X A , X B ] : X A , X B E g r } => g r , = {X 3 , X 4 } (3.14) 

and furthermore, it's second derived algebra reads: 

g r „ = {[X A ,X B ] : X A ,X B G g r ,} => g r „ = {0} (3.15) 

Thus, the group G r is solvable since the g r „ is zero. As it is evident A 3 ,A 4 ,y2 
generate an Abelian subgroup, and we can, therefore, bring them to their normal form 
simultaneously. The appropriate transformation of the dependent variables is: 



(3.11) 



(3.12) 



7n = e 

712 = e 

713 = e 
722 = e 



U1+2UQ 



UI+U2+U4+U6 



itl+2 114 



(3.16) 



733 = e Ul (1 + e 2 " 6 u 3 2 + 2 e U2+U4+m u 3 u 5 + e 2u4 u 5 2 ) 
In these coordinates the generators Y 2 , X A assume the form: 



Y = — 

2 dui 



X = 9 



V 9 



Sua dii3 



(3.17) 



Except of the parametrization ( 13.16)1 there is also another one achieving the same result 
(13.17)1 . which simply attributes a - sign to 712 and therefore any solution later described 
will remain valid under this change. 

Evidently, a first look at ()3.16)) gives the feeling that it would be hopeless even to write 
down the Einstein equation. However, the simple form of the generators (|3.17j) ensures 
us that these equations will be of first order in the functions in, it 3 and M5. 
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3.1 Description of the Solution Space 

Before we begin solving the Einstein equations, a few comments for the possible values 
of the functions Ui, % = 1, . . . , 6 will prove very useful. 
The determinant of 7 Qi g, is 

det[ laP ] = e 3«i+2(«4+« 6 ) ( X _ e 2« 2 j ( 318 ) 

so we must have u 2 < . 

The transformation from the 7' s to the it' s, becomes singular when 712 = 0, since 
the function u 2 equals to 

, n h hi(7n 722) , q in x 

u 2 = ln(|7i 2 |) . (3.19) 

So two cases are naturally arising, according to whether 712 is different or equal to zero. 
If 7i2 7^ the two linear constraint equations, written in the new variables (j3.16|) . give 

E t = -e U(i (e U(i u 3 + e U2+U4 u 5 ) = (3.20) 
E 2 = -- e UA ( e U2+ue u 3 + e U4 u 5 ) = (3.21) 

This system admits only the trivial solution, since the determinant of the 2x2 matrix 
formed by the coefficients of ii 3 ,u 5 becomes zero only for the forbidden value u 2 = 0. 
We thus have 

u 3 = fa, u 5 = k 5 (3.22) 

Now, these values of u 3 ,u 5 make 713,723 functionally dependent upon 711,712,722 (see 
(I3.16J1 ). It is thus possible to set these two components to zero by means of an appro- 
priate constant automorphism. 

In the case 712 = we can again bring simultaneously into normal form the correspond- 
ing X 3 , X4, Y 2 . The appropriate change of dependent variables is given by: 

(gUl+2u6 Q g«l+2«6 ^ \ 

e «i+2«5 g Ul -«4+«5 (3.23) 

e «l+2« 6% eUl -u 4 +u 5 e «! ^ _|_ e -2u 4 _|_ e 2u 6 U 2^J 

In these variables all three linear constraint equations can be integrated, yielding: 

£i = 0^ -e 2 " 6 ii 3 = => u 3 = k 3 (3.24) 

E 2 = => ~ e - U4+U5 (u 4 + m 5 ) = =>• u 5 = k 5 - u 4 (3.25) 

E 3 = 0^ -2 e 2u4+2u6 u 3 u 3 + u i + u 5 + 2e 2u4 u 6 = 0^u 6 = k 6 (3.26) 
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. Again, these values imply that a constant automorphism suffices to set the (13) and 
(23) components of the scale-factor matrix to zero, i.e. to put it into diagonal form. We 
have thus reached a first important conclusion, that is: 

Without loss of generality, we can start our investigation of the solution space for Type 
III vacuum Bianchi Cosmology from a block-diagonal form of the scale-factor matrix ( 
and, of course, zero shift) 

(7n 7i2 \ 
7i2 722 (3.27) 
733/ 

Note that this conclusion could have not been reached off mass-shell, due to the fact that 
the time-dependent Automorphism ()3.4j) does not contain the necessary two arbitrary 
functions of time in the (13) and (23) components ( besides the fact that all the freedom 
in arbitrary functions of time has been used to set the shift to zero). As we have earlier 
remarked, since the algebra ()3.13j) is solvable, the remaining (reduced) generators X\, X 2 
(corresponding to diagonal constant automorphisms) as well as Y 2 continue to define a 
Lie-Point symmetry of the reduced EFE's and can thus be used for further integration 
of this system of equations. 



3.1.1 Case I: 712 = 

The remaining (reduced) automorphism generators are 

d d 
X\ = 2711- — , X 2 = 2722- 



22 




The appropriate change of dependent variables which brings these generators -along with 
Y 2 - into normal form, is described by the following scale-factor matrix : 

laP = I e U2+u * I (3.28) 

In these variables the first two linear constraint equations are identically satisfied, while 
the third reads £3 = =>• — 2«x = =>• u\ — kx- Substituting this value of u\ into the 
quadratic constraint equation E we obtain the lapse function 

N 2 = ^-e Us u 3 (2u 2 + 3u 3 ) (3.29) 
lb 

. Now, substitution of u% = k\ and the above value for the lapse iV 2 into the spatial 
EFE's results in the single, independent equation : 

(ii 2 + M3) (2M3M2 - 2m 2 -u 3 + 2ulu 3 + 3-u 2 + 5-U2W3) (3.30) 

This equation is, as expected from the theory, of the first order in U2,u 3 . Notice that 
this result could have not been reached had we chosen any particular time gauge, such 
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as N 2 = F(u 2 ,u 3 ,t) : Not only u 2 ,u 3 ,t would appear in the Spatial EFE's, but also 
the number of independent such equations would have been increased to 2. This remark 
should not be taken as a negative view for complete gauge fixing, but rather as pointing to 
the fact that keeping the gauge freedom into the game helps manifesting the symmetries 
of the system and eventually solving the equations. 

Equation ()3.30|) is readily integrated, leading to two different space-times according to 
which parenthesis is set to zero. If the first is made to vanish, i.e. u 2 = k 2 — u 3 , the 
ensuing line-element is the known (Type III) cosmological disguise of Minkowski space- 
time ([E]): 

ds 2 = -^e U3 u 2 3 dt 2 + \e U3 dx 2 + e kl+u *~ 2x dy 2 + e k2 dz 2 (3.31) 
16 ' 4 

the constants being of course absorbable by the constant automorphisms and a shift in 

If the second parenthesis of ()3.30|) is set to zero, i.e. u 2 = k 3 — ^f- + ln(l + k 2 e 2 ), we 
obtain an equivalent form of the Type III member of the known Ellis-MacCallum family 
of solutions (0,115!): 

ds 2 = K 2 ( e^ii 2 df2 + £U3 ^ + eU3 _ 2x dy2 + ^ _ dz2 \ 

V 4(e^-l) J 

where again we have used constant automorphisms and a shift of u 3 to take outside of 
the metric an overall constant. We can be assured that the constant is essential either 
by checking that indeed the metric inside the parenthesis does not admit a homothetic 
Killing vector field or, more primarily, by finding an invariant relation between curvature 
and higher derivative curvature scalars which explicitly involves k. For metric ()3.32j) one 
such invariant relation is: 

18 Ql _ 2 n _ dKLMN d n — n u KLMN R 1"? ii\ 

K , Lj\ — K tiKLMN-, V2 — U -K -tiKLMN [6.66) 



(Q2- 



>1:AQ1;B yj 



where capital Latin letters denote space-time indices ranging in the interval (0-3), the 
semicolon stands for covariant differentiation, and the □ for the covariant D'Alebertian. 
This relation, being a constant scalar constructed out of the intrinsic geometry ( the 
Riemmann tensor and its covariant derivatives ), characterizes , along with many others 
that can be found, this metric: It will be valid for any equivalent, under general coor- 
dinate transformations, form of ()3.32j) . It is also noteworthy to observe that in both 
the above line-elements the arbitrary function of time appears; This is because the 
number of symmetry generators matches the number of scale-factors (both are 3), so 
that the system of spatial EFE's is reduced to first order without any choice of time. In 
the case of a block- diagonal scale-factor matrix one of the four scale factors will have 
to play the role of time before the corresponding system can be reduced. Lastly, met- 
ric ()3.32|) admits, except of ()3.3j) . a fourth Killing vector field acting on the surfaces of 
simultaneity, namely 

U = -2yd x + (e 2x - y 2 )d y (3.34) 
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There is thus a G 4 symmetry group acting (of course, multiply transitively) on each V 3 of 
this metric, with an algebra having the following table of (non- vanishing) commutators: 



[6, 6] = 6, [6,6] = -2 6, [6, 6] = 6 (3.35) 

However, it is interesting to note that we have not imposed the extra symmetry from 
the beginning , but rather it emerged as a result of the investigation process. 

3.1.2 Case II: 7l2 ^ 

The remaining (reduced) automorphism generators are 

d d ^ d d 
X 1 = 2711- h 7127^ — , a 2 = 7l2 - h 27 22 - 



. r\ ' I _1_Z O 7 Z / ±Z r\ 1 / ZZ r~\ 

O711 C?7l2 O712 <^722 

The appropriate change of dependent variables which brings these generators -along with 
Y 2 - into normal form, is now given by: 

(gUl+2«4 gMi+U2+«4 Q 

e «l+n 2 +n 4 e «i+2«2 M3 I (3.36) 

The generators are now reduced to 

Odd , 
Y 2 = —,X 2 = —,X 1 = — 3.37 

OU\ OU2 OUi 

indicating that the system will be of first order in the derivatives of these variables. 
The remaining variable u 3 will enter, (along with ii 3 , u 3 ) explicitly in the system and 
is therefore advisable (if not mandatory) to be used as the time parameter, i.e. to effect 
the change of time coordinate 

t — > u 3 (t) = s, ui(t) — > ui(t(s)), u 2 {t) — > u 2 (t(s)), u 4 (t) — > -u 4 (t(s)). (3.38) 

This choice of time will of course be valid only if u 3 is not a constant. We are thus led 
to consider two cases according to the constancy or non-constancy of this variable. 

The case u 3 = k 3 

The determinant of the scale-factor matrix becomes det[y a p\ = e 3n i+ 2 («2+«4) _|_ ^ 
We thus have k 3 > 1. The two linear constraint equations are identically satisfied, while 
the third yields 

171 ii 2 + (l-2k 3 )u 4 u 2 
E 3 = ^ - — = u 4 = k 4 + 



2(1 -k 3 ) 2k 3 
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Inserting these values of u 3 , u 4 into the quadratic constraint equation we obtain the 
following lapse 

= e"* (-I + /C3) (3(l-2£; 3 ) 2 n 2 + 8A; 3 (-l + 2fc 3 ) t^ 2 + 4A^ 2 ) 

4(l-2fc 3 ) 2 (-3 + 4A; 3 ) 

Use of these values of u 3 , 114, (N) 2 in the spatial EFE's results, as expected, in a system 
which is of the first order in the unknown variables Ui, ii2- The coefficient of 112 in 
£33 = is 

2e Ml faui ((-1 + 2fc 3 ) U1 + U2) 

3(1-2 k 3 ) 2 u 2 + 8k 3 (-1 + 2 k 3 ) mi ii 2 + 4 A; 3 u 2 . 

and can be safely regarded different from zero, since the possibilities iii = 0, u 2 = 
(1 — 2 k 3 ) u\ easily lead (through E 33 = itself) to zero and negative lapse, respectively. 
We can thus solve E 33 = for 112 and substitute into E12 = which becomes 

e k 4+Ul+ _ 1+2k3 ^ ^ i + 2il2 ) ((_ 3 + 6 fe 3 ) u l + 2 (3-2fc 3 ) tx 2 ) _ Q 

6 - 20 k 3 + 16 fc 3 2 

Again, the second parenthesis in the numerator leads to zero lapse, leaving us with the 
only alternative iii = —2iL2 => u\ = 2k\ — 2u2 which indeed satisfies all spatial EFE's. 
Finally, inserting these values of u 3 , w 4 , U\ in the lapse 1)3.39]) and the scale-factor matrix 
()3.36|) we obtain the following line-element (after using the constant automorphisms and 
a shift in U2 to purify the metric from the absorbable constants) : 

ds 2 = -A 2 d? + ^dx 2 + e~ 2x ^ x dy 2 + dz 2 + 2e~ x £ 2X dy dz (3.40) 

4 2A — 1 

where the constant A is related to k 3 by k 3 = ^=4 =>- < A < ~ and we have adopted 
the time gauge e~ U2 = £ for simplicity. 

This metric is an equivalent form of a solution originally given by Siklos ^Hj and repro- 
duced in [15] . An overall multiplicative constant has been omitted from ()3.40j) since it 
admits the following Homothetic Killing vector field (£# qab = A* 9ab) 

^4 +(i - 2A) 4 + 4 

It also admits three more Killing vector fields ( except ()3.3)1 ) acting on space-time, namely 
£ 4 = e ~fx d^ + ^e-ih d x 



___ . , ,2Ay 3L A(A — 1) 4x-i — 1 \- l 1 . 1 

£5 = e ^yd^ + — -e 2A ^ + ^ — r e 2X £ d v 



2A-1 



2Xz 
J 

-A(2A- l)e~^id z 



e 2 a z -\ — e 2\d x — Xe ^ x £ 2A+1 d y 
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The £5 field breaks down for A = |, in which case the valid expression is 

y 3 (-2x + io g (0) -yg 

The first of these is null £4 gf 5^ = and covariantly constant £f. B = 0, signaling 
that the metric is a pp-wave. Consequently all scalar curvatures, constructed by forming 
scalar contractions of tensor product of the Riemmann tensor and its covariant deriva- 
tives of any order (such as Qi, Q2 in ()3.33j) ). vanish identically (see e.g. [Ej). This raises 
the interesting question of how can we be certain that the constant A is essential. An 
answer can be found in terms of equalities between tensors -constructed out of the Riem- 
mann tensor and its covariant derivatives- that hold true in these space-times [T8].|19j. 
For metric ()3.4()j) such a relation is : 

R A B° D RaECF;G;H = . : 9 , „ > ~ R A B ° D;E RaFCG;H (3.42) 

— AA Z + 2A + 2 

By the quotient law the expression of A in the right-hand side of this relation is a 
scalar function, and being a constant it can not change value under any coordinate 
transformation; thus A can not be altered by such a transformation and is, therefore, 
essential. 

The algebra of the six killing fields, ()3.3|) . (j3.41j) has the following table of non 
vanishing commutators: 

Ei, 61 =6 = & = [6^5] = Es,&] = ^6(3.43) 

There is an isotropy group G2 of null rotations emanating from this algebra, which is 
easily seen by taking a linear combination of these fields: 



n = 6 ^ = 6 ^6 = 6 (3-44) 
e.g. [F 1 ,y 2 ]=y 2 and \Y l ,Y< i \ = Y z 

The space (being pp-wave) does not obviously have curvature singularities, it thus 
seems to be geodesically complete and is of Petrov Type N. 

The case w 3 7^ k 3 

The function w 3 is now a valid choice of time and det[) a p\ = e 3u i+ 2 ( M 2+M 6 ) ^_]_ _|_ s ^ 
implies the range (1, +00) for the new time s. The only non-vanishing linear constraint 
equation E 3 = yields 

w 4 = / - — ^— ds + (3.45) 
/ zs 1 
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while the quadratic constraint equation E = gives the lapse 



(Nf 



4(1 -2 s) (-3 + 4 a) 



[2 (2 s - l) 2 mi + 3 (2 s - l) 2 (s - 1) m 2 



(3.46) 



+ (4 s - 2) ti 2 + 8 s (s - 1) (2 s - 1) iii ii 2 + 4 s (s - 1) m 2 ,] 



If we insert these values (iV) 2 , u 4 into the spatial EFE's they become the following 
polynomial system of first order in ui, ii 2 



ii\ = (l Mi u\ iii) 



/1\ 

ill 



, u 2 = (luiiijul) A 2 



/1\ 

ii 2 

ill 



(3.47) 



/ 



Ai 



i 



4s 2 -7s+3 

2s-3 
4s-3 

6 s (s-1) 
4s-3 



Is 



4s 2 -7s+3 



16s 2 (s-1) 
'8s 2 -10s+3 



Ss 2 -10s+3 



is (2 a-3) (s-1) n 
8s 2 +10s-3 U 



0/ 



(3.48) 



A, 



/ 



-4s+2 
4s 2 -7s+3 



-6s+3 
4s-3 



-8s+5 24s 2 -50s+18 8 s (2 s-3)(s-l) \ 

8s 3 -18s 2 +13s-3 8s 2 -10s+3 8s 2 -10s+3 ' 



12 s 

-2s+3 

6 s (s-1) 
4s-3 





16 s 2 (s-1) 
8s 2 +10s-3 













/ 



(3.49) 



Due to the form of Ai, A 2 (their components are rational functions of the time s), 
system (13.47(1 can be partially integrated with the help of the following Lie-Baklund 
transformation 



ui(s) 



(2 s - 3) tanr(s) - 2 s (8 s 2 - 10 s + 3) r(s) 
4 s y/s^T (4 s - 3) 



2s - 1 



u 2 (s 



8s(4s-3) y/(s - l) 3 



(2 (-4 s + 3) - 1 + 3 (s - 1) tanr(s) 



(3.50) 



+2s{s- 1) (4s-3)f(s)) 
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resulting in the single, second order ODE for the variable r(s) 



\/s=T\ 2 (-16s + 6)v / s^I+(5s-3) tanr 

tanr r H = r 

2 J 2s(4s-3) y/(s-l) 

-9(s-l) 2 tan 2 r + 18(s-l) 3 / 2 tanr + 4 s (4 s - 3) 
+ 8 s 2 (4s-3) 2 (s- l) 3 /2 t 3 ' 5 ) 

This equation contains all the information concerning the unknown part of the solution 
space of the Type III vacuum Cosmology. Unfortunately, it does not posses any Lie- 
point symmetries that can be used to reduce its order and ultimately solve it. However, 
its form can be substantially simplified through the use of new dependent and indepen- 
dent variable (p,u(p)) according to r(s) = ± arcsin -^== , s = 3 ^ p p Z^ > P > § thereby 

obtaining the equation 

u = ± . u 2 = w „ n r 3.52 

v/(6p-10) (p 2 -u 2 -l) (6p-10)(p 2 -« 2 -l) V ; 

with the corresponding lapse 

W 2 - 8 (3p - 5)^-^-1) ^ < 3 - 53 > 

(u = j 4 ) and the scale-factor matrix is given by ()3.36|) after insertion of ([3.45)1 . % = 

s = an d the transformations of ui, 112 that led to u. Independently of the way we 

have reached this result, one can check (through an algebraic computing facility such as 
Mathematica) that the line element thus described is indeed a solution of all the EFE's, 
provided of course ()3.52|) is satisfied. One can also check that it does not admit any 
Homothetic or null, covariantly constant vector field. Therefore, the two independent 
constants of the general solution to 1)3.52)1 along with a multiplicative constant will com- 
prise the expected three essential constants of the general Type III vacuum Cosmology: 
The general algorithm for calculating this number when a space time gauge has been 
chosen (usually zero-shift and unit lapse), in which case the constraints must be viewed 
as restrictions on the initial data, reads as : 

12 ( for the six components of jap) -1 ( for the time reparameterization covariance )- 
number of independent constraints -dimension of Automorphism Group. 
When a space-time gauge has not been fixed, i.e. when constraints are being viewed as 
symmetry generators, the relevant counting is given by 



D = 2 x (number ofjap) 

— 2 x (number of linear constraints) 

— 2 x (Quadratic constraint) 

— (number of parameters of outer-Aut) 
-(n) 
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where n = dim(inner-Aut) - number of independent linear constraints. 

In our case the number of independent linear constraints is 3, and the dimension 
of the inner- Aut is 2, so n = — 1. The constants that appear at the outer- Aut are 2 
and obviously the number of 7 Q/ g is 6. Thus, the expected maximal number of essential 
constants is indeed 3, by both ways of counting . 

Despite the relatively simple form of ()3.52|) . its general solution is, to the best of our 
knowledge, not known. However, we have managed to obtain a partial solution in the 
parametric form 

4(l + 2e 2 ^) 3 / 2 
3(l + e 2 «) 2 

which makes the functions u±, U2, read as 

w i(0 — ki + £ + In cosh£ 

u 2{0 = k 2 + In sech^ — - In (cosh2^ + 2) 

W4(£) = &4 — In coshf; + — In (cos/i2£ + 2) 

and the lapse (iV) 2 = 6 jr^if^rp^ - The ensuing metric, after the usual purification with 
the constant automorphisms and a shift in £, is given by : 

ds 2 = k 2 - - \ c ( d£ 2 + — coshi dx 2 + e- 2x+ Hcosh2£ + 2) sech£ dy 2 
\ 4(2e 2 ? + l) 4 v ' 

+ e ? sechi dz 2 + 2e~ x+i sech£ dy dz^j (3.54) 

As we have already remarked, this metric does not admit a Homothety and therefore 
the constant k is essential. It does not satisfy the invariant relation ()3.33|) . and it is not 
a pp-wave. Therefore we conclude that it is inequivalent to ()3.32|) or ([3.40)1 . This mono- 
parametric family belongs to the Kinnersley vacuum solutions jO] . It is quite interesting 
that it also admits a fourth killing vector field 

£ 4 = -lQyd x + (e 2x -8y 2 )d y -2e x d z (3.55) 

which produces with (|3.3J) the following table of (non-vanishing) commutators: 

[6,6] = £1, = -16^, = & (3.56) 

The isotropy group inferred from the above algebra (see the last commutator) is a 
G\ spatial rotation. 

Curiously enough, this algebra is equivalent to ()3.35|) as a simple scaling of £1, 6 by 
2^/2 shows. Of course, the multiply transitive character of the action of the underlying 
group on the corresponding V3 's allows for these, and thus for the space-times in which 
they are embedded, to be inequivalent. 

Again, the extra symmetry emerged in the course of the investigation of the solution 
space. Of course it must have something to do with the particular nature of the solution, 
but it was not set starting point. 
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3.2 Preview for other Bianchi Types 

The method described in the previous sections can be applied to other Types as well. 
The general pattern is s similar to that of Type III : The pp-wave solutions ( for Types 
admitting such geometries) occupy one part of the solution space, the other known 
solutions reside on another part, and the unknown part of the solution space is always 
described by an ODE strikingly similar to (J3.52J1 . namely : 

2_ (-1+U 2 ) 2 

Details will be included in a forthcoming work. As indicative examples we give the form 
of the ODE for Types IV and VII h : 
Type IV 

C_l + ^2)2 

u 2 = x v - ' k = —6, A = 6 (3.58) 



(k + X P ) (p 2 -u 2 -l) 



Type VII t 



-1 + ii 2 ) 2 



■,2 T" " J „ _ a , 4 



and of course 
Type III 



4 Discussion 

When one is trying to solve Einstein's Equations in cosmology, one has to deal with 
a nonlinear system of coupled, ordinary differential equations. The strategy that is 
frequently used, is to simplify the system by choosing, a convincing form for the scale 
factor matrix, usually obtained by an a priori assumption of extra symmetry (e.g. 7^ = 
diag(a(t),b(t),c(i)) and then try to solve it, hoping to find some solution. Clearly, 
this procedure can by no means guarantee access to the full space of solutions for the 
problem at hand. In this work we have presented a method for solving Einstein's Field 
Equations in the case of vacuum Bianchi Geometries. The main idea is to consider the 
Group of constant Automorphisms, which emerges as the residual freedom left after the 
time dependent A .I.D.'s f!2.8|) . (j2.9|) have been used to set the shift N a to zero, as a 
Lie point-symmetry of the EFE's. In a step-by-step procedure one can bring some of 
the generators of this group in normal form and simplify the rest, thereby reducing the 
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order of the system of equations. Which of the generators, and how, can be utilized 
in each step depends upon the characteristics of their Lie Algebra (abelian, solvable 
etc.). It is also important that the information gained at a particular level must be used 
and, in fact, may be vital for the implementation of the next step. The method is, by 
construction, sweeping out all possible solutions, since no ad-hoc assumption has been 
made. Therefore, if successfully applied to a given Bianchi Cosmology, it will result in 
the cartography of the entire space of solutions. 

The successful application of the procedure to Bianchi Type III resulted in the recov- 
ery of all known solutions without prior assumption of any extra symmetry ( (|H.H2j) . (|H.4()|) ). 
the enclosure of the entire unknown part of the solution space into a single, second order 
ODE in terms of one dependent variable ([3.52)1 . and a partial solution to this ODE. It 
is of interest that the solution space is naturally partitioned into three distinct discon- 
nected pieces. Of great importance may be considered the fact that a strikingly similar 
ODE describes the unknown part of the solution space for other lower Bianchi Types. 
For Types VIII, IX there remain no rigid automorphisms after the shift has been set to 
zero and the constant rotations have been used to diagonalize the scale-factor matrix. 
However there is the scaling symmetry Y 2 that can serve as a starting point. This issue, 
along with the presentation of the detailed cartography for the lower Types is in our 
immediate scopes. Finally the method can be extended towards either the inclusion of 
matter content, or in 4+1 Spatially Homogeneous Cosmologies. 
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